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Abstract 

This paper addresses the performance of space-time coding over fading channels with impulsive 
noise which is known to accurately capture network interference. We use the symmetric alpha stable 
noise distribution and adopt two models which assume dependent and independent noise components 
across receive antennas. We derive pairwise error probability (PEP) of orthogonal space-time block 



C/3 



' codes (STBC) with a benchmark genie-aided receiver (GAR), or the minimum distance receiver (MDR) 



m 



X 
J3 



which is optimal in the Gaussian case. For general space-time codes we propose a maximum-likelihood 
(ML) receiver, and its approximation at high signal-to-noise ratio (SNR). The resulting asymptotically 



(N 

' optimal receiver (AOR) does not depend on noise parameters and is computationally simple. Monte- 



Carlo simulations are used to supplement our analytical results and compare the performance of the 



receivers. 
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I. Introduction 

The additive Gaussian noise model has long been used because it produces simple and 
tractable mathematical models which are useful for gaining insight into the underlying behavior 
of communication systems. As the physical reality of most practical channels demonstrate much 
more sophisticated effects such as bursts and impulses, which arise as a consequence of man- 
made activity such as automobile spark plugs [1], microwave ovens [2], and network interference 
[3]-[8], the Gaussian noise model may not be accurate. Such environments are also observed 
in urban and indoor channels as well as underwater acoustic channels [9], [10]. Therefore, 
impulsive noise which captures these physical effects should be considered. In such wireless 
environments, the performance is degraded both by fading and impulsive noise. To combat 
fading, antenna arrays are often used, giving rise to multi- input multi-output (MIMO) systems. 
Space-time coding has been used as one of the powerful diversity techniques in MIMO systems. 

A number of performance analyses of STBC have been reported in the literature where the 
noise is Gaussian (see e.g., [11]-[13]). Recently, some works in the area of STBC in the presence 
of impulsive noise have also been reported. Performance of space-time diversity/coding for power 
line channels with Middleton Class-A noise model was studied by simulations in [14]. In [15] the 
code design criteria and the PEP upper bound were derived over a fading channel with Middleton 
Class-A noise. Subsequent work in [16] provided a closed-form expression for symbol error rate 
(SER) of orthogonal STBC (OSTBC) when the noise follows a Gaussian mixture model. 

Symmetric ce-Stable (SaS) distributions are an important class of noise distributions which 
can successfully model a number of impulsive noise processes. Studies [3]-[8] show that, in a 
multi-user network with power-law path loss, the multiple access interference results in a SaS 
distribution, when the interfering nodes are scattered according to a spatial Poisson point process 
(PPP). In [17], the performance evaluation of a MIMO system in SaS noise was performed by 
simulation with no closed-form expression for the error probability. Subsequent works in [8] and 
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[18] provided closed form expressions for the bit error rate (BER) of linear diversity combining 
schemes for SaS noise environments in single-input multi-output (SIMO) environments. In [19], 
[20], the optimal linear receivers for SaS noise were studied in SIMO systems. To the best of 
our knowledge there is no analysis of MIMO systems over fading channels with SaS noise. To 
close this gap in the literature, our goal is to design receivers for, and analyze the effect of SaS 
noise on space-time coded systems. While the receivers derived herein apply to all space-time 
codes, the (PEP-based) performance analysis holds for OSTBCs. 

Throughout this paper, we use (•)^ for Hermitian, ( j^ for transpose, diag(x) for a diagonal 
matrix with elements of x along the diagonal, || || for the Frobenius norm for matrices and 
Euclidean norm for vectors, Ai(-) for the i^'^ largest eigenvalue of a matrix, to denote the 
real part, ^f-} to denote the imaginary part. Also, we use Ea,b{C) to denote the expected value 
of the random variable C with respect to the distributions of the random variables A, B. Finally, 
we write f{x) — 0{g{x)) as x — >■ a to indicate that limsup^..^^ \ f{x)/g{x)\ < oo. 

II. System Model 

We consider a wireless communication system where the transmitter is equipped with Nt 
antennas and the receiver with N^. antennas. We consider the following standard MIMO flat- 
fading channel model: 

Y = VPHS + W (1) 

where Y is the N^xTg received signal matrix, and Tg is the length of the transmitted data block; 
H is an A^rX^t matrix, with independent and identical distributed (i.i.d.) circularly symmetric 
complex Gaussian entries with mean zero and variance 1 ; the average transmitted power at each 
transmitting antenna is denoted by the scalar p; S is the NfXTg transmitted data block, which is 
transmitted from a codeword set S with equal probability; W is the N^xTg additive impulsive 
noise matrix, with elements that have a SaS distribution, as explained next. 
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We first introduce real valued SaS random variables, which will later be used to define its 
complex counterpart used in this paper. A real valued (not necessarily symmetric) a-stable 
random variable, w Sa{a, jj) has a characteristic function given by [21], [22] 
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a e (0,2] is the characteristic exponent, /3 e [—1, 1] is the skew, cr e (0, oo) is the scale and 
// e (— oo, oo) is the shift parameter. When /3 = 0, w has a symmetric distribution about When 
/9 = and /i. = 0, w is a SaS random variable. When a = 2 and /9 = 0, w is Gaussian, which 
is the only SaS random variable with finite variance. Since the Gaussian case is widely studied, 
we focus on a e (0, 2) throughout. When a —\ and = 0, w is said to be standardized [23, 
pp. 20]. Any SaS random variable w Sa {cr, 0, 0) can be written as compound Gaussian, i.e., 
of the form w = VAG, where A and G are independent, with A ~ Sa/2 ([cos(7ra/4)]'/'^,l,0) 
is positive skewed a-stable random variable and G ~ ^2 (ex, 0, 0) is Gaussian random variable 
with mean zero and variance 2cr^ [21, pp. 38], [23, pp. 20]. 

Although a closed-form expression for the PDF of SaS random variables exists only for a few 
special cases (e.g. Gaussian (a = 2) and Cauchy (a = 1)), asymptotic expansions for a e (0, 2) 
are well known as w — > oo: 

Uw) = a(l + /3)CaW-"-^ + 0(w-2"-i) (5) 
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where the constant := r(Q;) sin(7rQ;/2)/7r [22]. Additionally, if w ~ Sa{o;/3,0), the com- 
plementary cumulative distribution function (CCDF) of w satisfies the asymptotic relation as 
A — >■ oo: 

P(^/; > A) = Ca(7"(l + /3)A-" + 0(A-'"). (6) 

In the following, we will briefly introduce two noise models (Model I and II) which assume 
dependent and independent noise components across antennas. In both Model I and II, the Tg 
columns of W, Wi, Wr^, in (1) are independent. 

• Under Model I, we assume := W2,fc, ^Af,.,*;]^ is a complex isotropic SaS 
random vector, defined as 

= v^(G^ + iG^,) (7) 

where the scalar random variable ~ 3^/2 is independent 

of G§ and G;^ which are Gaussian random vectors with i.i.d. elements which have 
mean zero and variance cr^. This is a good assumption when the receiving antennas 
are influenced by the same physical process creating the impulse, thereby making the 
Ak of each branch the same. This might, for example, be an accurate model for a 
multi-antenna system where the antenna elements spaced closely. Mathematically, it is 
not difficult to see that in this case u'i fc,u'2,fc, .. , wat^ ^ will be statistically dependent 
[23, pp. 83]. 

• Under Model II, the j, k element of W is given by 

[Wkfe = v^(Gf,+iG},,) (8) 

where Aj^^: Gfk and G^ f. are distributed as in Model I, but are i.i.d., and [W]j_fe is the 
(j, k) element of matrix W. 
In both Model I and II, wj^k has a unity scale parameter (a — 1), since any scale is subsumed 
in p in (1). It can be shown that only the moments of order a or less exist for any SaS random 
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variable [21, pp. 22], as a result of which the conventional definition of SNR holds only for the 
Gaussian case (a = 2). However, with a slight abuse of terminology, we will refer to p as the 
SNR, even when a < 2, since p quantifies the relative scale of the signal versus the noise. 



We assume throughout that the channel H is known at the receiver. Under Model I, we start 
with the GAR for which are assumed known at the receiver at each time k — 1, ...,Ts. The 
GAR is optimal when Ak are known, so that its performance can serve as a benchmark for any 
practical receiver that does not have this knowledge. 

A. Genie-aided Receiver 

The GAR maximizes the posterior probability and hence minimizes the probability of error, 
when H and Ai,..., At^ are known. In the following, we are going to derive the decoding rule. 
To express in matrix form, we define A — diag([l/\/]47, I/a/AtJ). Right multiplying (1) 
by A, we obtain: 



so that the product WA has i.i.d. CA/^(0, 1) entries. Since the elements of WA are now white 
Gaussian and the codewords are equally likely, the optimal decision rule is to minimize the 
Euclidean distance: 



To express the PEP that S is transmitted and S' is received for the GAR in (10), we follow 
the derivation in the Gaussian noise case and obtain. 



III. Receiver Design and Performance 



YA = ^HSA + WA 



(9) 



S = argmin||YA - y/pHSAf. 



(10) 



s 




(11) 



Using (11) and Craig's representation of the Q function. 




(12) 
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(11) can be expressed as follows: 



P(S^ S'|H,A) = - 



TT 



/p||H(S-S')Aii2 
exp ' 



Taking expectation with respect to H and A, we get 



Nt 
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de. 



Nr 



(13) 



(14) 



where B := (S — S')AA^(S — S')^. Using (14), we can show that the code design criterion 
under SaS noise remains the same as the Gaussian noise case as follows. To obtain the maximum 
diversity order, we need B to be a full rank matrix for any realization of A in (14). Since A 
is diagonal with nonzero diagonal elements, it is a full rank matrix. Therefore, if the codeword 
different matrix S — S' is full rank, B is guaranteed to be a full rank matrix. 

When S — S' is square and unitary which is satisfied by e.g., the Alamouti code [24], the 
eigenvalues satisfy Ai(B) = 1/Ai. Substituting in (14), using the statistical independence of Ai, 
and taking expectation, we show in Appendix A that, as p — > oo 



P(S^S') 



1 r(^) 



2 



r (i + f)r(iv,.^f) 
r (1 - f ) T(N,) 



p 



=:GGAK{Nt,Nr,a) 



2 



+0 {p 



(15) 



Using {Gc ■ p)'^"^ expression to present PEP, we can define the diversity order, Gd, and the 
coding gain, Gc, from the PEP. The coding gain is defined as the amount that bit energy or 
signal-to-noise power ratio can be reduced under the coding technique for a given bit error rate. 
In (15), the Gc is GcARi^t, and the G^ is aNt/2. The implications of (15) are interesting, 

because it suggests that the diversity order depends on the number of transmit antennas, Nt, and 
the noise parameter, a. However, the number of receive antennas, Nr, does not contribute to the 
diversity order. This is due to the fact that the noise is not i.i.d. across antennas in Model I. 
In order to investigate the behavior of the coding gain as a function of Nr, by differentiating 
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the natural logarithm of the coding gain with respect to Nr, we get 

^ logGoAR = -- [Nr - f ) - V' (Nr)] (16) 

where ^{x) :— is the digamma function as defined in [25, pp. 258-259]. In (16), since 

V'(x) is a monotonically increasing function for x > 0, the term inside the brackets is negative 
Va e (0,2). Therefore, the coding gain is a monotonically increasing function of N^.. So, even 
though Nr does not contribute to diversity, it does improve the coding gain. Regarding the 
analysis of GcAKiNt, Nr,a) in (15) with respect to Nt, it is shown in Appendix B that the 
coding gain is a monotonically decreasing and convex function of Nf. 

For Model II the GAR can also be derived by using the Hadamard product with A which is 
a matrix with ( j, k) element 1 / ■\fA~k- However, its performance is not tractable. 

B. Minimum Distance Receiver 

The MDR, which is optimal over Gaussian noise minimizes the Euclidean distance: 

S = argmin||Y - ypHS||l (17) 
s 

Note that unlike the GAR in (10), the MDR does not depend on A. We now derive the PEP for 
the MDR. Define E := H(S - S')/||H(S - S')||, and let e^-fc be the (j, A;) element of E. The 
PEP and its upper bound for the MDR are given by: 




pIIH(S-S') 


p 










P||H(S-S 


)ll 






(18) 
(19) 

(20) 
(21) 
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where A^nax '■— max^^lfc is the maximum value among Ai, ...,At^. In (20) we used the fact 
that ||E|| = 1, and in (21) we used (12). Taking expectation with respect to H and Amax, the 
following upper bound on the average PEP is obtained: 



Eh,a„.P (S ^ S'\A, 



TT Jo 



Nt 



Nr 



Si 



p Ai(C) 



(22) 



where C := (S - S')(S - S') . When S - S' is square and unitary, we can rewrite (22) as 
follows: 



£h,4„.,P(S-»S'|A 



TT Jo 
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4 sin^ d Amo 

Taking expectation with respect to A^ax, we show in Appendix C that, as p — > oo 



P(S^ S') < 



N, r {^) r {NrN, - f ) 
20F r (1 - f ) r(Ar^7vo ' 



(24) 



=:GMDR(Art,Ar^,Q) 

Equation (24) suggests that the diversity order is always a/2 regardless the number of antennas 
which is reduced compared to the GAR where it was aNt/2. 

The behavior of the coding gain as a function of can be obtained from the derivative given 

by 



d 



logG 



MDR 



2N, 



a 



a 



(25) 



In (25), since ^jJ{x) is a monotonically increasing function for x > 0, we can verify the term inside 
the brackets is negative Va e (0,2). Therefore, the coding gain is a monotonically increasing 
function of Nr. Next, by differentiating the log-coding gain with respect to Nt, we get 



— — log Gmdr = 

oNt a 



— + Nr 

Nt 



lP (^NrNt - I) - V' {NrNt) 



(26) 



It can be shown numerically that the GuDRiNt, Nr, a) monotonically decreases with Nt when 
a e (0, ceo) for some constant ckq. Unlike the GAR, in case of the MDR the number of transmit 
antennas, Nt, does not contribute to the diversity order. Hence when a G (0, ao) the performance 
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of MDR will be worse as Nt increases. Intuitively, the reason for the deterioration in performance 
is that when a is small, the sum of independent noise samples do not "average out" like it does 
when the noise has a finite variance. In other words, when a is small enough the performance 
bound of MDR suffers from increased transmit antennas! On the other hand, the coding gain is a 
monotonically increasing function of Nt when o; e (cti, 2) for some constant cti. In other words, 
when a e (ai, 2) the coding gain increases as the number of transmit antennas increase. When 
a e (tto, cti), the coding gain is a concave function of Nt. The values of cto and cti depend on 
Nr (e.g. when Nr = 1, ao ~ 1.333 and ai 1.799). 

For Model n the PEP of MDR can be derived by using Amax '■= maxj ^ Aj^k in (19). Following 
the same derivation, the PEP of MDR for Model II is obtained by multiplying GuDRiNt, Nr, a) 
in (24) with Nr^^°' which implies less coding gain and the same diversity order. This is in 
contrast with the GAR which will be shown in the simulations to have better performance under 
Model n compared to Model I. In conclusion, for SaS noise environments the conventional 
MDR receiver has poor performance especially for small a. 

C. Maximum Likelihood Receiver 

We introduce the optimal ML receiver for Model I and II. Firstly, the optimal ML receiver 
for Model I is given by 



fc=i 

where /q;(||x||) is a probability density function of amplitude distribution of d-dimensional 
multivariate isotropic stable random variables and is given by [26]: 



S = argmaxTT/„ (Hy^ - ^/pHsk\\) 



(27) 



argmax log (| | y ^ - ^/pUsk \ \ ) 



(28) 




(29) 



where r = llx 



^^Xf + \- and J,y{-) is the Bessel function of order u. 
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In case of Model II, complex symmetric a-stable random variables are independent in both 
space and time. Thus, we can modify the optimal ML receiver for Model II as follows: 

Ts Nr 

S = argmax^ J]log/,(||[Y],-fc - ^[US],^). (30) 
^ k=i j=i 

Since faif) cannot be expressed in terms of closed-form elementary functions, these ML re- 
ceivers are seen to be computationally complex, and dependent on the noise parameters a and 
a. We now consider receivers that perform nearly optimally, with the advantage of reduced com- 
plexity and not requiring knowledge of noise parameters, when compared to the ML receivers. 

D. Asymptotically Optimal Receiver 

To simplify (27), we use the expression for the tail of /a(-) in [26] 

= ^,.E^n(a + 2)^^a + i)m^^,,,., ^ 0(.-<-«)) (31) 
as r ^ oo, where we note that a2'^ ^^"(W^) r((a+2)/2)r((«+d)/2) q ^ ^ . dominant 

-Tra/Z 1 (a/ 2) ' ^ 

term of (31) in (27) and simplifying, we get 

Ts 

S = argminTT||yfc - ^Hsfcll (32) 

fe=i 

Ts 

= argmin V'logllyfe - V^Hsfell. (33) 

Using same approach as Model I, we can modify the asymptotically optimal receiver for Model 
II as follows: 

Ts Nr 

S = argmin^^log||[Y],,fc - ^^^[HS],, (34) 
^ k=i j=i 

The resulting receivers are asymptotically optimal at high SNR and relatively simple. 

A few comments about complexity of the ML receiver and AOR follow. In (28) and (33), we 
need to evaluate matrix norms. The only difference between (28) and (33) is that the equation 
(28) needs to evaluate the metric in (29) additionally. In (29), it is needed to evaluate an 
elementary function, a special function (i.e., the Bessel function) and an integration of these 
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functions for each candidate codeword S. Instead of evaluation of (29), we can alternatively 
use a lookup table for the numerical values of (29). Such a lookup table would have sizable 
memory requirements since a lookup table would contain values for each of the a. and a values 
corresponding to the noise parameters. For example, if the sizes of quantized a. and a values are 
and A^^ respectively, we need the N^. ■ N^r entries in the table. In addition to these kinds of 
high computational complexity, the ML receiver also requires to estimate a and a values of SaS 
noise. However in case of the AOR which performs within a tenth of a dB of the ML receiver 
which will be shown in Section IV, we do not need to evaluate the equation (29) and estimate 
the a and a values. 

Therefore we propose to use the AOR for impulsive noise due to its relatively low complexity 
and its reasonable performance. Though our analysis is based on the receivers for Model I, it 
will be similar in case of Model 11. We note that the asymptotically optimal receivers in (33) 
and (34) are additive and therefore can be used in conjunction with the Viterbi algorithm when 
S is a codeword on a trellis. 

IV. Simulations 

In this section, we verify our results through Monte Carlo simulations. In our simulations, 
we assume that a = 1.43, which corresponds to the value estimated in [27] for modeling radio 
frequency interference in laptop receivers. We also consider a "highly impulsive" scenario, with 
a = 0.5, which corresponds to a path loss exponent of 2 /a = 4 in an environment where the 
interfering nodes are scattered according to a PPP on a two-dimensional plane [8]. 

A. Performance Results under Model I 

We show in Fig. 1 the performance bound of GAR for Alamouti code with Nt = 2, Nr = 1 
over highly impulsive noise with BPSK. We calculate the BER union bound using the PEP of 
GAR in (15). We also plot the upper bounds for the MDR obtained using (24). In Fig. 1, we 
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also show the simulated BER results of Alamouti code for GAR, MDR, ML receiver, and AOR. 
Comparing between theoretical and simulated results, we observe the diversity orders of GAR 
and MDR are aNt/2 and a/2. We also observe the performance gap between ML receiver and 
GAR is about 1.3 dB at 10~^ BER. We also found the performance for AOR which does not 
need the noise parameters shows a difference less than a tenth of a dB to the ML receiver. 

In Fig. 2, we show the performance of Alamouti code with A^^ = 2, A^^ = 2. It is noted that 
the diversity orders do not change even though the number of receiver antennas increases in 
accordance with our theoretical result. In this case, the ML receiver and AOR are seen to be 
within 0.6 dB of the GAR. 

In the following, we show the performance of Alamouti code over impulsive noise with 
a — 1.43. In Fig. 3, we show the theoretical and simulated BER with Nt — 2,Nr — 1. 
The performances with Nt = 2, Nr = 2 are shown in Fig. 4. Under the less impulsive noise 
environment with a — 1.43, we observe that the diversity orders of GAR and MDR are also 
aNt/2 and a/2 which are in line with our theoretical results. It is also observed the performances 
for ML receiver and AOR are within 2.5 dB of the GAR at 10~^ BER, as suggested by Fig. 3. 
In Fig. 4, the ML receiver and AOR are seen to be within 1 dB of the GAR. 

B. Performance Comparison between Model I and 11 

In Fig. 5, we compare the simulated performances of Alamouti code over highly impulsive 
noise under Model I and II. Under Model II, we can observe that the diversity order of GAR 
will be larger than that of Model I, because additional diversity can be obtained due to the 
independence of the noise in the space domain. However, the diversity order of MDR does not 
change even under Model II. We can also observe the performance difference for AOR and ML 
receiver is less than a tenth of a dB. Additionally, we show the simulated performances over 
impulsive noise with a = 1.43 under Model I and II in Fig. 6 where we observe the diversity 
order of GAR of Model II is larger than that of Model I and the diversity orders of MDR are 
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always ajl. It is also observed the performance difference for AOR and ML receiver is less 
than a tenth of a dB. 



V. Conclusions 

In this paper, we considered a SaS noise model for MIMO fading channels, and discussed 
different receivers. In SaS noise environments, the diversity order depends on the noise param- 
eter, a, and noise correlation model. Under Model I, we derived the diversity order for the GAR 
and MDR. The maximum possible diversity order of GAR is shown to be a benchmark for 
any receiver, given by aNt/l. The MDR, though simple, is vulnerable to impulsive noise: the 
diversity order is always a/2 regardless the number of antennas. Under Model II we have seen 
that the diversity order for GAR will be larger than that of Model I. In contrast, for MDR the 
diversity order is q;/2 also for Model II. 

Since the GAR is impractical to implement, we are motivated to use the ML receiver. However, 
the ML receiver is computationally complex and requires knowledge of the noise parameters. 
Thus, we also develop an asymptotically optimal receiver, which performs near optimally at 
high SNRs and does not require the noise parameters. Since the conventional MDR has poor 
performance, the usage of the MDR should be avoided in SctS noise environments. 



Appendix A 

Since {A}i=i,...,Art are i.i.d., each term of the product in (14) has the same expected value, 
which results in 
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where A represents any of the random variables Ai. To simplify the expectation in the RHS of 
(35), recall that A ~ ,Sa/2([cos(7rQ;/4)]2/", 1, 0), so that 
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(36) 
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The PDF /a/2(^) as suggested by (5) as — > oo is given by 



(37) 



Substituting (37) in (36), we get 



\AAsm'^6 + p 

ar(f)r(iv.-|) 



C„/2Cos(^)^-(^+t)+o(^-(^+")) 



2 r(i- f)r(iv,) J \Asm'^e 
r(i + f)r(iv,-|) ^ ^ p 



r (1 - f) r{Nr) J Usin^e 

Plugging (38) in (35), and using the binomial expansion, we get 



(38) 



P(S^ S') = - 
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r(i + f)r(iv,-f)\ /p 



r (1 - f ) T{Nr) 

Solving the integral in (39), (15) follows. 



(9 



Nt 
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Appendix B 

To prove the coding gain, GoARiNt, -/V^, a), is a monotonically decreasing and convex function 
with respect to Nt, we will show a stronger statement which states that the coding gain is a 
logarithmically completely monotonic (cm.) function which means that the derivatives of the 
logarithm satisfy: 



(-ir 



d 

m 



logGGAR{Nt,Nr,a)>0 



(40) 



for n e Z+. Letting aNt/2 = x in the coding gain of (15), it suffices to show that 

r(x + i) = 

is a logarithmically cm. function. Taking logarithm in (41), we get 



(41) 



fix) = logh(x) = -log2x/7r+ -logr(a; + 1) - -logT (x + -] . (42) 
X X X \ 2 J 
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Using Leibnitz' rule, [u{x)v{x)]^''^ = J2k=o {l)u^''Kk)v^''~''Kk), in each of the last two tenns 
in (42), we obtain 
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X 



k=0 



n\ f 1 



k \k 



(n-k) 



[iogr(x + i)]^')- 



logr (x + 



X 



g{x) 



(43) 



where 



g{x) :=\og2^+\ogV{x+l)-\ognx+-)+Y.^^^^l^^'-'\^^^^ lx+- 

k=l ' k=l ' ^ 



and '0(")(a;) is the polygamma function as defined follows [25, pp. 260]: 
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1 - e-* 
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Jo 



dt. 



(44) 



(45) 



The proof will be complete when we show g{x) > for x > since that would make (43) 
positive. The first derivative of g{x) can be expressed as follows: 



g'ix) = 



[-l)''x' 



ni 



iP^^'^x + 1) - iP^""^ (x+^ 



in) 



Using (45), we conclude 



x"- n\ jQ 



re-^* dt > 0, 



(46) 



(47) 



1 - e-* 

since (^ ^ ^^_~-t * ^ > for i > 0. Thus, the function g{x) is increasing and g{x) > g{0) > on 
(0, oo), which implies (x) > and n = 0, 1, 2, . . . . Thus, h{x) is a logarithmically 

cm. function. Since a logarithmically cm. function is also cm. [28], h{x) is a cm. function, 
which in turn has convex and decreasing functions as a special case. Therefore, it is proved the 
coding gain is a monotonically decreasing and convex function with respect to Nt. 

Appendix C 

Using (6) and order statistics, we can find the PDF of Amax as follows: 

= |r,2C./.cos (^) (l - 2C„,, cos (^) .-f + O (.-t««)-) . 

(48) 
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If > 2 and is equal to Nt, plugging (48) in (23) and using binomial expansion, we can get 
as follows: 



Ax sin^ 9 
yAxsin^O + p 

Nt~l 



NrNt 



aNt_ 
2 



Nt-l . , \ 

k=0 J 



E 

k=G 



(-1)' 



a(l+fc) 



r ( NrNt - 



(r (i-f))^'+'^r(iv,7v, 



p 



a{l + k) 



+ O [p —) 



2 r (l - f ) T{NrNt) y V4sin2 9 
' NtV{l + f^V{NrNt-f) 



r (1 - I) r(iv,iVi) 

Plugging (49) in (23), we get 



V4sin 9/ 



(49) 



P (S ^ S') 



^U^^(^^f)^v^•^'-^)\(^^^e^oiP)-^. (50) 



ttJo V r (1 - f ) r(7v,7v,) 



Solving the integral in (50), (24) follows. 
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Fig. 1. Performance comparison of GAR, MDR, ML receiver, and AOR over a channel with highly impulsive noise (a = 0.5) 
with iVt = 2 and iV,. = 1 




Fig. 2. Performance comparison of GAR, MDR, ML receiver, and AOR over a channel with highly impulsive noise (a = 0.5) 
with Nt =2 and Nr = 2 



January 12, 2013 



DRAFT 



21 



10" 



10 



10"' 



10"' 



tr 

LU 

m 



10 



10" 



10"' 



10 



(a,Nt,Nr)=(1 .43,2,1 ) Alamouti scheme 




-* — GAR (simulation) 
— IVIDR (simuiation) 
iVIL (simulation) 
□ AOR (simulation) 

GAR (upper bound) 

IVIDR (upper bound) 



10 



15 



20 



25 
SNR 



30 



35 



40 



Fig. 3. Performance comparison of GAR, MDR, ML receiver, and AOR over a channel with moderately impulsive noise 
(a = 1.43) with iV* = 2 and iV^ = 1 
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Fig. 4. Performance comparison of GAR, MDR, ML receiver, and AOR over a channel with moderately impulsive noise 
(q = 1.43) with Nt=2 and Nr = 2 
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(a,Nt,Nr)=(0.5,2,2) Alamouti scheme 
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Fig. 5. Performance comparison of GAR, MDR, ML receiver, and AOR over a channel with highly impulsive noise (a = 0.5) 
with iVt = 2 and iV,. = 2 under Model I and II 
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Fig. 6. Performance comparison of GAR, MDR, ML receiver, and AOR over a channel with moderately impulsive noise 
(q = 1.43) with Nt=2 and Nr = 2 under Model I and II 
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